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j—{ . Abstract 

> : 

. We review the nuclear physics input necessary for the study of the collapse of 

f-*) I massive stars as precursor to supernova explosions. Recent theoretical advances for 

' the calculation of the relevant weak-interaction processes and their influence in the 

■ collapse and postbounce dynamics are discussed. Future improvements are expected 
! to come with the advent of radioactive ion-beam facilities. 

S ■ 1 Introduction 

-t— > ■ 

. 

I Stars with masses exceeding roughly 10 Mq reach a moment in their evolution 

^ ' when their iron core provides no further source of nuclear energy generation. 

■ At this time, they collapse and, if not too massive, bounce and explode in 
spectacular events known as type II or Ib/c supernovae. These explosions 
mark the formation of a neutron star or black hole at the end of the life of the 
star and play a preeminent role in the nucleosynthesis and chemical evolution 
of the galaxy. The evolution in the core is determined by the competition of 
gravity, that produces the collapse of the core, and the weak interaction, that 
determines the rate at which electrons are captured and the rate at which 
neutrinos are trapped during the collapse. 

The early phases, known as presupernova evolution, follow the late-stage stel- 
lar evolution until core densities just below 10^° g cm~^ and temperatures 
between 5 and 10 GK are reached. Stellar evolution until this time requires 
the consideration of an extensive nuclear network, but is simplified by the 
fact that neutrinos need only be treated as a sink of energy and lepton num- 
ber. This is no longer valid at later stages of the collapse: as the weak in- 
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teraction rates increase with the increasing density, the neutrino mean free 
paths become shorter so that the neutrinos eventually proceed through all 
phases of free streaming, diffusion, and trapping. An adequate handling of the 
transitions between these transport regimes necessitates a detailed time- and 
space-dependent bookkeeping of the neutrino distributions in the core. During 
collapse, only the electron neutrinos, z/e, arc important because the positron 
abundance, which would lead to electron antincutrino emission by capture 
on neutrons, is very low under electron-degenerate conditions. Later in the 
evolution the electron degeneracy is partially lifted and in addition to the 
electron flavor neutrinos also heavy neutrinos, v^, 1/^, and their antiparticles, 
are usually included in numerical simulations of core collapse and postbounce 
evolution. 

Advantageously, the temperature during the collapse and explosion are high 
enough that the matter composition is given by nuclear statistical equilibrium 
without the need of reaction networks for the strong and electromagnetic 
interaction. The transition from a rather complex global nuclear network, 
involving many neutron, proton and a fusion reactions and their inverses, 
to a quasi-statistical equilibrium, in which reactions within mini-cycles are 
fast enough to bring constrained regions of the nuclear chart into equilibrium, 
to finally global nuclear statistical equilibrium is extensively discussed by [1]. 

The crucial weak processes during the collapse and postbounce evolution 
are [2,3,4,5,6] 



p + e ^n + Ue (1) 

n + + Ue (2) 

{A, Z)+e-^ {A, Z-l) + u, (3) 

{A,Z)+e+^{A,Z+l) + u, (4) 

jy + N^iy + N (5) 

N + N^N + N + U + P (6) 

u + (A,Z)^u + (A,Z) (7) 

z/ + ^ z/ + (8) 

u+{A,Z)^iy + {A,Zy (9) 

e+ + e^^z/ + z/ (10) 

{A,zy^{A,Z) + iy + iy (11) 

Ve^e^^lJ.,T^H,T (12) 



Here, a nucleus is symbolized by its mass number A and charge Z, N denotes 

cither a neutron or a proton and v represents any neutrino or antincutrino. 
We note that, according to the generally accepted collapse picture [7], elastic 
scattering of neutrinos on nuclei (7) is mainly responsible for the trapping, 
as it determines the diffusion time scale of the outwards streaming neutrinos. 
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Shortly after trapping, the neutrinos are thermahzed by energy downscatter- 
ing, experienced mainly in inelastic scattering off electrons (8). The relevant 
cross sections for these processes are collected in appendix C of Ref. [2] (for 
an update see Ref. [6] in this volume). Reactions (1) and (3) are equally 
important, as they control the neutronization of the matter and, in a large 
portion, also the star's energy losses. Due to their strong phase space sensitiv- 
ity (~ E^), the electron capture rates increase rapidly during the collapse as 
the density (the electron chemical potential scales like ~ p^^^) and the tem- 
perature increase. In the postbounce phase, it is neutrino-nucleon scattering, 
Eq. (5), that provides the main opacity, and lepton capture on nucelons, eqs. 
(1) and (2), that are responsible for the dominant creation and absorption of 
electron flavor neutrinos. While reactions (4), (9) and (11) have not yet been 
routinely included in dynamical models, nucleon-nucleon bremsstrahlung [8,9], 
Eq. (6), and UgP^ annihilation [10,11,12], Eq. (12), have replaced the traditional 
electron-positron annihilation, Eq. (10), in the role of the dominant source for 
II and T neutrinos. 

Numerical simulations of core collapse have a long history (see Ref. [7] and ref- 
erences therein). After the observation of the closeby supernova 1987A in the 
Magellanic Cloud, simulations in spherical symmetry reached a climax around 
1990. At that time it was established that the neutrino transport scheme has 
to be performed in separate energy groups to account for the strong (oc E^) 
energy dependence of the neutrino cross sections and it was realized that 
the thermalization of the neutrino distribution function by neutrino-electron 
scattering plays an important role in the deleptonization and core dynamics 
[2,13,14,15,16,17]. Most simulations relied on the flux-limited diffusion equa- 
tion for the neutrinos. The accuracy of this approximation in the collapse 
phase was evaluated by comparison to approaches where the full Boltzmann 
transport equation was solved [3]. It took another decade until these more 
sophisticated schemes also mastered the postbounce phase in general rela- 
tivistic space-time [18,4]. The narrowing technical uncertainties in spherically 
symmetric simulations and the accessible information about the detailed neu- 
trino distribution functions sparked new interest to improve the nuclear and 
weak interaction physics ingredients. These efforts coincided with the devel- 
opment of large-scale shell model calculations (see Ref. [19] for a review) that 
allowed for a more sophisticated and accurate treatment of electron capture 
rates for the presupernova and collapse phases [5,20,21]. Moreover, the shell- 
model calculations could be validated [22] against the experimental data on 
GT distribution measured on charge-exchange experiments. The advent of new 
radioactive ion beam facilities and high resolution experimental techniques 
based on the (i, ^He) [23,24] and (d, ^He) [25,26] charge-exchange reactions 
will contribute to an improved treatment of the electron capture rates used in 
simulations. 

No spherically symmetric simulation of core collapse and postbounce evolu- 
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tion with accurate neutrino transport produced an explosion [27,28,29,30,31], 
neither by the prompt bounce-shock mechanism, nor by the delayed neutrino- 
heating mechanism. Despite the mentioned progress in simulation technique 
and input physics, it became very evident that the complicated dynamics of 
the delayed supernova explosion can no longer be expected to be captured in 
spherically symmetric simulations. It has been shown that the heating region 
is unstable to convective overturn which significantly improves the efficiency 
of the neutrino heating behind the stalled bounce-shock [32,33,34]. However, 
later investigations with two-dimensional hydrodynamics and more accurate 
(particularly energy-dependent) neutrino transport still failed to reproduce 
vigorous supernova explosions [35,36]. Furthermore, suggestions that an en- 
hanced neutrino luminosity due to neutron finger convection in the protoneu- 
tron star [17] could increase the heating behind the shock are also not con- 
firmed by recent analytical analysis [37]. Hence, more ingredients to a realistic 
supernova model must be missing or inadequately be accounted for. There are 
however indications, that an additional twist might suffice to at least produce 
an explosion, like for example a very asymmetric deformation of the standing 
accretion shock [38,39] or effects from magnetic fields [40,41,42,43]. 

In this review we will discuss nuclear physics ingredients to which core col- 
lapse supernovae are sensitive. The manuscript is organized as follows: first we 
discuss the role of weak interactions in the presupernova evolution. Then we 
comment on the main hydrodynamical features of core collapse based on a re- 
alistic nuclear equation of state. We will see that the dynamics is governed by 
entropy changes and deleptonization, which we then investigate in detail with 
respect to traditional and new nuclear physics ingredients. Finally, we discuss 
the infiuence of the recently improved collapse physics on the post-bounce 
evolution. 



2 Presupernova evolution 

The main weak interaction processes during the final evolution of a massive 
star are electron capture and beta decays. Its determination requires the cal- 
culation of Fermi and Gamow-Teller (GT) transitions. While the treatment of 
Fermi transitions (important only for beta decays) is straightforward, a correct 
description of the GT transitions is a difficult problem in nuclear structure. 
In the astrophysical environment nuclei are fully ionized, so one has contin- 
uum electron capture from the degenerate electron plasma. The energies of 
the electrons are high enough to induce transitions to the Gamow-Teller reso- 
nance. Shortly after the discovery of these collective excitation Bethe et al. [44] 
recognized its importance for stellar electron capture. This process is mainly 
sensitive to the location, fragmentation and total strength of the Gamow- 
Teller resonance. The presence of a degenerate electron gas blocks the phase 
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space for the produced electron in beta decay. Then, the decay rate of a given 
nuclear state is greatly reduced or even completely blocked at high densities. 
However, due to the finite temperature, excited states in the decaying nucleus 
can be thermally populated. Some of these states are connected by large GT 
transitions to low-lying states in the daughter nucleus, which with increased 
phase space can significantly contribute to the stellar beta decay rates. The 
importance of these states in the parent nucleus for the beta decay was first 
recognized by Fuller, Fowler and Newman (FFN) [45,46,47,48], who coined 
the term "backresonances" (see figure 1). 

Electron capture 




(Z+1 ,A) 

Fig. 1. The figure shows schematically the electron-capture and beta-decay processes 
in the stellar environment. Electron capture proceeds by Gamow- Teller transitions 
to the GT^ resonance. In the case of beta decay both the Fermi and Gamow- Teller 
resonances are typically outside the Qp window and hence not populated in the beta 
decay. Due to the finite temperature in the stellar environment, excited states in the 
decaying nucleus can be thermally populated. Some of these states have large GT 
transitions to low lying states in the daughter nucleus. These states in the decaying 
nucleus are called "backresonances" 

Over the years many calculations of weak interaction rates for astrophysical 
applications have become available [49,50,51,52,53,54,55]. For approximately 
15 years, though, the standard in the field were the tabulations of Fuller, 
Fowler and Newman [45,46,47,48]. These authors calculated rates for electron 
capture, positron capture, beta decay and positron emission plus the associ- 
ated neutrino losses for all the astrophysical relevant nuclei ranging in mass 
number from 21 to 60. Their calculations were based upon an examination of 
all available experimental information in the mid 1980s for individual transi- 
tions between ground states and low-lying excited states in the nuclei of in- 
terest. Recognizing that this only saturated a small part of the Gamow- Teller 
distribution, they added the collective strength via a single-state represen- 
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tation whose position and strength was parametrized using an independent 
particle model. 
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Fig. 2. Comparison of the shell-model GT-|- distribution (lower panel) for ^^V 
with the high resolution (d, ■^He) data [56]. The shell-model distribution includes 
a quenching factor of (0.74)^. 

Recent experimental data on GT distributions in iron group nuclei [57,58,59,60,61,62,63,64], 
measured in charge exchange reactions [65,66], show that the GT strength is 
strongly quenched (reduced), compared with the independent-particle- model 
value, and fragmented over many states in the daughter nucleus. Both effects 
are caused by the residual interaction among the valence nucleons. An accu- 
rate description of these correlations is essential for a reliable evaluation of 
the stellar weak-interaction rates due to the large dependence of the available 
phase-space on the electron energy, particularly for the stellar electron-capture 
rates [45,20]. The shell- model is the only known tool to reliably describe GT 
distributions in nuclei [67]. Indeed, Ref. [22] demonstrated that the shell- model 
reproduces very well all measured GT_|_ distributions (in this direction a pro- 
ton is converted to a neutron, as in electron capture) for nuclei in the iron 
mass range and gives a very reasonable account of the experimentally known 
GT_ distributions (in this direction a neutron is converted to a proton, as in 
P decay). However, the limited experimental resolution (~ 1 MeV) achieved 
by the pioneering (n,p)-type charge-exchange experiments did not allow for a 
detailed determination of the fragmentation of the GT strength in individual 
states. Very recently, high-resolution GT+ distributions measured at KVI, via 
the {d, ^He) reaction, have become available for two iron group nuclei, ^^V [56] 
and ^^Ni [68]. The experimental data for ^^V are compared in figure 2 with a 
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shell-model calculation using the KB3G interaction [69] . 



Several years ago, it was pointed out that the interacting shell model is the 
method of choice for the calculation of stellar weak-interaction rates [55,70,71,72,73] . 
Following the work of Ref. [67], shell-model rates for all the relevant weak 
processes for srf-shcll nuclei {A = 17-39) were calculated in Ref. [74]. This 
work was then extended to heavier nuclei {A = 45-65) based on shell-model 
calculations in the complete p/-shell [20,75]. Following the spirit of FFN, the 
shell model results have been replaced by experimental data (energy positions, 
transition strengths) wherever available. 



Ref. [20] compares the shell-model based rates with the ones computed by 
FFN. The shell-model rates are nearly always smaller than the FFN ones at 
the relevant temperatures and densities. The differences are caused by a reduc- 
tion of the Gamow- Teller strength (quenching) compared to the independent- 
particle-model value and a systematic misplacement of the Gamow- Teller cen- 
troid (mean energy value of the Gamow- Teller distribution) in nuclei depend- 
ing on the pairing structure. In some cases, experimental data that were not 
available to Fuller, Fowler and Newman, but could be used now, led to signif- 
icant changes. 
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Fig. 3. Comparison of the center values of (left), the iron core sizes (middle) 
and the central entropy (right) for 11-40 Mq stars between the WW models and 
the ones using the shell model weak interaction rates (LMP) [76]. The lower parts 
define the changes in the 3 quantities between the LMP and WW models. 

The influence of the shell-model rates in the late-stage evolution of mas- 
sive stars has been investigated in Ref. [76,77] repeating the calculations of 
Ref. [78], keeping the stellar physics as close to the original studies as possi- 
ble. The new calculations incorporated the shell-model based weak-interaction 
rates for nuclei with mass numbers A = 45-65, supplemented by the srf-shcll 
nuclei rates from Ref. [74]. The earlier calculations of Woosley and Weaver 
(WW) used the FFN rates for electron capture and an older set of beta decay 
rates [51,52]. Figure 3 illustrates the consequences of the shell model weak 
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interaction rates for presupernova models in terms of the three decisive quan- 
tities: the central electron-to-baryon ratio Y^, the entropy, and the iron core 
mass. The central values of Y^. at the onset of core collapse increased by 0.01- 
0.015 for the new rates. This is a significant effect. For example, a change 
from Ye = 0.43 in the WW model for a 20 Mq star to Ye = 0.445 in the new 
models increases the respective Chandrasekhar mass by about 0.075 Mq. We 
note that the new models also result in lower core entropies for stars with 
M < 20 Mq, while for M > 20 Mq, the new models actually have a slightly 
larger entropy. The iron core masses are generally smaller in the new models 
where the effect is larger for more massive stars (M > 20 Mq), while for the 
most common supernovae (M < 20 Mq) the reduction is by about 0.05 Mq. 
[We define the iron core as the mass interior to the point where the composi- 
tion becomes at least 50% of iron group elements {A > 48)]. This reduction of 
the iron core mass appears to be counterintuitive at first glance with respect 
to the slower electron capture rates in the new models. It is, however, related 
to changes in the entropy profile during silicon shell burning which reduces 
the growth of the iron core just prior to collapse [76]. 




Time till core collapse (s) Time till core collapse (s) 

Fig. 4. Evolution of the Ye value in the center of a 15 Mq star (left panel) and a 
25 Mq star (right panel) as a function of time until bounce. The dashed line shows 
the evolution in the Woosley and Weaver models (WW) [78], while the solid line 
shows the results using the shell-model based weak-interaction rates of Langanke 
and Martmez-Pinedo (LMP). The most important nuclei in the determination of 
the electron-capture rate for the calculations adopting the shell model rates are 
indicated at different times. 

To understand the origin of these differences it is illustrative to investigate 
the role of the weak-interaction rates in greater detail. The evolution of Ye 
during the presupernova phase is plotted in figure 4. Weak processes become 
particularly important in reducing Ye below 0.5 after oxygen depletion (~ 10'' 
s and 10^ s before core collapse for the 15 Mq and 25 Mq stars, respectively) 
and Ye begins a decline which becomes precipitous during silicon burning. 
Initially electron capture occurs much more rapidly than beta decay. As the 
shell model rates are generally smaller than the FFN electron capture rates, the 
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initial reduction of is smaller in the new models; the temperature in these 
models is correspondingly larger as less energy is radiated away by neutrino 
emission. 

An important feature of the new models is shown in the left panel of figure 4. 
For times between 10^ and 10^ s before core collapse, Ye increases due to the 
fact that beta decay becomes temporarily competitive with electron capture 
after silicon depletion in the core and during silicon shell burning. This had 
been foreseen in Ref. [79]. The presence of an important beta decay contri- 
bution has two effects. Obviously it counteracts the reduction of Yg in the 
core, but equally important, beta decays are an additional neutrino source 
and thus they add to the cooling of the core and a reduction in entropy. This 
cooling can be quite efficient as often the average neutrino energy in the in- 
volved beta decays is larger than for the competing electron captures. As a 
consequence the new models have significantly lower core temperatures than 
the WW models after silicon burning. At later stages of the collapse, beta 
decay becomes unimportant again as an increased electron chemical potential 
drastically reduces the phase space. 

We note that the shell model weak interaction rates predict the presuper- 
nova evolution to proceed along a temperature-density- trajectory where 
the weak processes are dominated by nuclei rather close to stability. Thus it 
will be possible, after next generation radioactive ion-beam facilities become 
operational, to further constrain the shell model calculations by measuring 
relevant GT distributions for unstable nuclei by charge-exchange reaction, 
where we emphasize that the GT+ distribution is also crucial for stellar f3- 
decays [71]. Figure 4 identifies those nuclei which dominate (defined by the 
product of abundance times rate) the electron capture during various stages 
of the final evolution of 15 Mq and 25 Mq stars. An exhaustive list of the 
most important nuclei for both electron capture and beta decay during the 
final stages of stellar evolution for stars of different masses is given in Ref. [77] 

In total, the weak interaction processes shift the matter composition to smaller 
Yg values (see Fig. 4) and hence more neutron-rich nuclei, subsequently affect- 
ing the nucleosynthesis. Its importance for the elemental abundance distri- 
bution, however, strongly depends on the location of the mass cut in the 
supernova explosion. It is currently assumed that the remnant will have a 
larger baryonic mass than the iron core, but smaller than the mass enclosed 
by the oxygen shell [80]. As the reduction of Y^ occurs mainly during sili- 
con burning, it is essential to determine how much of this material will be 
ejected. Moreover, as show in refs. [81,82,83], weak interactions determine the 
composition of the innermost supernova ejecta. Another important issue is 
the possible long-term mixing of material during the explosion [84]. Changes 
of the elemental abundances due to the improved weak- interaction rates are 
rather small because the differences with respect to FFN occur in regions of 
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the star which are probably not ejected. 



3 Core-Collapse and Bounce 

3. 1 Hydrodynamics 

At the onset of collapse of the progenitor model, the core dynamics becomes 
relevant. The dynamics of core collapse is an interesting example where nu- 
merical simulation stimulated the analytical understanding. Early spherically 
symmetric simulations (at that time adiabatic or based on leakage schemes) 
[85,86,87,88] suggested the separation of the collapsing material into an inner 
and an outer core. The inner core fails only marginally to be pressure sup- 
ported and the sound speed stays faster than the fluid velocity at all times. 
Under these conditions, the velocity develops to a linear function of the radius. 
However, because the matter density and sound speed decrease with increasing 
radius, there is a sonic point where the infall velocity exceeds the sound speed 
so that information about the collapsing inner core ceases to reach the outer 
layers. The outer layers therefore rather assume a free fall velocity profile in 
the rarefaction of the homologously collapsing inner core. In the following we 
review different factors that determine the size of the inner core and try to 
draw the connection to the initial energy imparted to the bounce-shock. 

Figure 5 shows typical scales for energy densities as a function of the rest 
mass density, p. The data has been taken from model G15 [31] which has 
been evolved with general relativistic Boltzmann neutrino transport. A time 
slice at bounce is shown. The order of magnitude of the energy density is 
given by the gravitational binding energy of mass shells (thick solid line). It 
is balanced by the kinetic energy (dash-dotted line) and the internal energy. 
The latter determines the fluid pressure, composed from the electron pressure 
(thin solid line), baryonic pressure (dotted line) and neutrino pressure (dashed 
line). The dominant pressure of the relativistic electrons scales with p"^/^ over 
a large density range (in Fig. 5 actually somewhat shallower due to the radial 
variations in the electron fraction, Y^). 

The essentially polytropic behavior of the equation of state allows an analytic 
investigation [90,91] of collapse. It shows that the mass of the inner core, Mic, 
is well approximated by 

Mie~(«:/«o)'^'Mo, (13) 

where M refers to the core mass and n to the coefficient in the polytropic 
equation of state p — Kp^ with 7 = 4/3. Values with index belong to the 
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Fig. 5. Energy densities as functions of rest 
mass density at bounce in model G15 [31]. 
Sliown are the gravitational energy den- 
sity, Sg, and the kinetic energy density, Ek- 
The total pressure splits into the electron 
pressure, Pe, the baryonic pressure, ph, and 
the neutrino pressure, p^. The latter has 
been determined by taking the second an- 
gular moment of the radiation intensity. 
Note that the dynamics barely depends on 
the neutrino pressure because of the domi- 
nant electron pressure at high density and 
the negligible coupling between neutrinos 
and matter at low density [89]. 
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Fig. 6. Velocity profiles at bounce (solid 
lines) . The time after maximum central 
density, tpb, is negative before bounce. 
The maximum infall velocity at 0.4 ms 
before bounce marks the edge of the 
inner core at that time. A circle is 
drawn at the sonic point. The later pro- 
files demonstrate how a pressure wave 
is launched at the center which runs 
through the inner core to turn into 
a shock front close to its edge. The 
dashed lines represent the sound speed 
as given by the Lattimer-Swesty equa- 
tion of state. 



marginally stable stage immediately before collapse. With the coefficient for 
the degenerate ultra-relativistic electron gas [92], 

one readily reproduces the conclusion that the mass of the inner core evolves 

proportionally to Y^. The evolution of the electron fraction is therefore the 
most significant ingredient to determine the mass of the inner core at bounce. 
Note that the pressure of trapped neutrinos does not reach the level of the 
electron pressure because their chemical potential is lower in accordance to the 
difference between the neutron and proton chemical potentials (cf. Fig. 7). 

The analytical approach becomes more complex for a general adiabatic index 
7 < 4/3. Nevertheless, it is possible to express the selfsimilar solution in 
terms of a time-independent density profile D{X), where X is a self-similarity 
variable. The inner core mass then scales with time, t, as [91] 
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M(X) = 47r f dxx^Dix). 
Jo 

G is the gravitational constant, M(X) is the mass enclosed inside the similarity 
variable X, and tbounce is the time at bounce. The analysis demonstrates that 
there is a dependence of the hydrodynamics of core collapse on the adiabatic 
index of the equation of state at subnuclear densities. Due to the dominance 
of the electron pressure, however, deviations from 7 = 4/3 cannot be large. 

The analytical analysis is only valid in the nonrelativistic limit. General rela- 
tivistic numerical simulations lead to a significantly smaller mass of the inner 
core [93]. Equation (15) predicts that the mass of the inner core decreases with 
time if the adiabatic index falls below the critical value for a stable core in 
Newtonian gravitation, i.e. 7 = 4/3. In general relativity, the effective grav- 
itational potential is deeper and the critical adiabatic index for stability is 
larger [94]. One may therefore argue that in the general relativistic case even 
the 7 = 4/3 from the ultra-relativistic electron gas falls below the critical 
adiabatic index. Hence, the inner core shrinks during collapse as suggested by 
Eq. (15). A rigorous argument, however, would have to rely on an analytical 
analysis of core collapse based on a general relativistic similarity variable [95]. 
Recent numerical simulations of core collapse specify a reduction of about 20% 
in the inner core mass due to general relativistic effects [29,18,96]. The size of 
the inner core is most important for the energetics of the core-bounce: 

The baryonic pressure displayed in Fig. 5 (dotted line) increases steeply at 
high density. Shortly above 10^^ g/cm^ the thermal pressure component is 
exceeded by the neutron degeneracy pressure. The plateau accross the phase 
transition from isolated nuclei to bulk nuclear matter at 10^^ g/cm^ is caused 
by the attractive nuclear forces. The phase transition with its many possible 
geometrical structures and statistical fragmentation is very interesting from a 
nuclear physics point of view [97,98,99,100]. However, a significant impact of 
the microscopic phenomena in this regime on the collapse and/or supernova 
dynamics remains to be demonstrated. At and above saturation density, re- 
pulsive nuclear forces dominate the stiffness of the equation of state. When 
the center of the collapsing core reaches these densities, collapse is halted by 
an outgoing pressure wave. This pressure wave travels through the inner core, 
where infall velocities are subsonic, and turns into a shock wave close to its 
edge. The matter in the inner core experiences a rather adiabatic pressure wave 
and remains at a low entropy ~ 1.4 kB per baryon. The shock wave in the 
outer core, however, heats matter to entropies larger than ~ 6 kB per baryon 
so that heavy nuclei are dissociated. If the bounce-shock were to dynamically 
propagate through the core to expel outer layers in a prompt explosion, it 
would have to provide the energy to dissociate the material between the edges 
of the inner core and the iron core at a rate of 1.5 x 10^^ erg per 0.1 Mq of 
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dissociated material (the shock additionally suffers from unavoidable neutrino 
losses) . 

The size of the inner core does not only set the mass that will be dissociated 
in the postbounce evolution, additionally it determines the initial energy of 
the bounce-shock. The rebound of a larger mass in the gravitational poten- 
tial imphes a larger initial energy imparted to the outgoing shock wave. This 
initial shock energy has been found to depend sensitively on the stiffness of 
the equation of state around nuclear densities [94,101,102,15,103,104,39]. Un- 
fortunately, most studies used the explosion energy of prompt hydrodynamic 
explosions to measure this effect. Their results are difficult to compare be- 
cause the final explosion energy depends on the physics applied to propagate 
the initial shock through the iron core. Some apparent contradictions in the 
conclusions of above-cited references suggest a few retrospective comments: 

The homologous core stays very large in adiabatic calculations when electron 
captures are neglected. Hence, energetic prompt explosions are obtained. The 
initial shock energy increases with a decrease of the adiabatic index in a poly- 
tropic equation of state [94]. A similar effect was found in a recent adiabatic 
comparison between two realistic equations of state with different incompress- 
ibilities at saturation density [104]. The period of the lowest oscillation mode of 
the inner core in the regime of the stiffened equation of state depends strongly 
on the difference of its mass from the maximum stable mass. If the mass of 
the inner core is much smaller than the maximum stable mass, the oscillation 
period is very small with respect to the sound crossing time in the outer part 
of the inner core. Infalling material piles up around a "hard" center until the 
outgoing compression wave reaches the edge of the inner core, where it turns 
into a relatively weak shock (see Fig. 6). If the maximum stable mass is close 
to the inner core mass, the oscillation period is larger and the oscillation is 
more efficiently excited. A larger mass bounces on a "soft" core and a more 
energetic shock is launched at its edge. As soon as electron captures during col- 
lapse are included, e.g. by a simple leakage scheme, the mass of the inner core 
shrinks and prompt explosions are more difficult to obtain. Only extremely 
soft equations of state lead to explosions in that case [101]. But sufficiently 
soft equations of state, producing a maximum stable mass close to the mass of 
the inner core, are not compatible with the significantly larger observed neu- 
tron star masses, unless a very extended phase transition provides softness at 
nuclear densities in combination with a stiff adiabatic index in the far supranu- 
clear regime [102]. Later investigations of the effect of the equation of state 
at bounce were consistent with the described trend, but ceased to produce 
energetic prompt explosions due to improvements in the neutrino transport, 
especially after the inclusion of neutrino electron scattering [13,14,15]. Few 
years later, the BCK equation of state [105] was replaced by the Lattimer- 
Swesty equation of state [89] in "standard" simulations of core collapse. Its 
maximum mass for the hot protoneutron star is around 2.4 Mq, much larger 
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than the mass of the inner core ~ 0.5 M0. The consequence is a very low ini- 
tial shock energy of 1.9 x 10^^ erg (measured in simulation G15 [31] according 
to the definition used in Ref. [15]) instead of the 5-10 x 10^^ erg obtained in 
previous studies with the BCK equation of state [15]. The dynamical shock 
stalls already at 3 ms after bounce. The ensuing expansion of the accretion 
front does not depend on the exhausted shock energy anymore. As we will 
briefly discuss in section 4, the expansion is determined by the accretion rate 
in the gravitational potential and the deleptonization rate of the accumulated 
matter. This explains why a comparison of different choices of incompressibil- 
ities at saturation density in Ref. [103] showed no discernible consequences for 
the shock propagation on a time scale of 50 ms. The study did not resolve 
the few milliseconds after bounce where transient effects of the compressibil- 
ity variations may manifest in the energy of the short-lived bounce-shock. An 
update on the comparison of three different equations of state on the longer 
time scale of postbounce simulations with accurate neutrino transport is given 
in Ref. [39]. 

In summary, the crucial element in core collapse calculations is the determi- 
nation of the mass of the inner core because it determines the point of shock 
formation at bounce and the initial energy that is imparted to the shock. The 
following paragraphs focus on the weak interactions which directly or indi- 
rectly determine the core mass through the evolution of electron fraction and 
entropy. 



3.2 Deleptonization 

The deleptonization and entropy changes during collapse are determined by 
the interplay between three different physical processes: (i) the transition of 
protons to neutrons by electron capture and neutrino emission (for example 
by reactions (1) or (3)), (ii) isoenergetic neutrino scattering (e.g. by reactions 
(5) or (7)), and (iii) neutrino thermalization (e.g. by reactions (8) or (9)). 
Each of these processes is rather straightforward to describe in a gas of free 
nucleons. In reality, however, the by far most abundant nuclear species in the 
cold collapsing matter are neutron- rich heavy nuclei (see Fig. 8) . While the de- 
scribed weak interactions basically stay the same, they rather involve nucleons 
bound in nuclei than free ones. This leads to the interesting expectation that 
the nuclear structure is probed in many nuclei that are difficult or impossible 
to explore under terrestrial conditions. In the following we first describe the 
principle importance of each of above processes for the deleptonization of the 
core before we get to their detailed quantitative description and dependence 
on nuclear structure in section 3.3. 

As the density, p, increases, the electron chemical potential increases with 
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fie oc p^/^ [7], see also Fig. 7 (thin solid line). The electrons fill higher energy 
levels and electron captures on free or bound protons become more and more 
likely. As long as the density is lower than 5 x 10^° g/cm^, the neutrinos escape 
freely and the deleptonization rate is determined by the electron capture rate. 
Figure 9 compares the electron capture reaction rates on free protons with 
the electron capture reaction rates on nuclei. The former is by an order of 
magnitude larger than the latter. But the nuclei are by orders of magnitude 
more abundant targets. It is difficult to decide from first principles whether 
electron captures on free protons or nuclei would dominate in core collapse 
(see for example [44,106,107]): Because the free proton fraction is very sen- 
sitive to small variations in the value, it can be huge if the Y^, value is 
somewhat larger than a reference value or it can be negligible if the Y^. value is 
somewhat smaller than a reference value. This had interesting consequences in 
core collapse simulations. As described above, progenitor models for different 
main sequence masses or with different sophistication of the included weak 
interaction physics may run into collapse with quite different electron fraction 
profiles. Collapse simulations with Boltzmann neutrino transport, however, 
showed that most of the initial variations in the inner core structure con- 
verged to the same trajectories during core collapse, leading to very similar 
positions of shock formation and nearly identical neutrino bursts [108,109]. 
These simulations used the standard treatment of electron capures [2] that 
inhibited the dominant electron capture on nuclei as soon as nuclei exceeded 
neutron number = 40. This leads to a more shallow deleptonization slope 
at a density ~ 2 x 10^° g/cm^ as illustrated in the left panel of Fig. 10. But 
the deleptonization is not slowed down for long because a high electron frac- 
tion at increasing density (where nuclei become more neutron rich) implies a 
steeply rising free proton fraction until the deleptonization resumes on a sim- 
ilarly effective time scale by captures on free protons. The Y^ value at which 
this happens is determined by the intrinsic properties of the nuclei at a given 
density and temperature and therefore independent from the initial electron 
fraction or proton abundance. Hence, the deleptonization of different mod- 
els converges toward a "norm" electron fraction trajectory. A convergence of 
this kind will always take place when the product of electron capture reaction 
rate and target abundance dominates for a nucleus whose abundance strongly 
depends on the electron fraction. 

This appears somewhat contradictory to results of earlier extensive studies of 
the effect of the free proton fraction on core collapse. Significant differences 
in the deleptonization and inner core size at bounce have been reported [14]. 
This result becomes understandable if one considers that the investigation has 
not been performed with different initial free proton mass fractions. Instead, 
different free proton mass fractions have been set by variations of the asym- 
metry energy in the equation of state. A smaller asymmetry energy lowers the 
electron fraction at which the free proton fraction becomes significant enough 
to perform electron captures effectively on the infall time scale. It therefore 
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leads to a lower "norm" trajectory for the electron fraction and to a smaller 
inner core at boimcc. However, if the asymmetry energy is varied together 
with the surface energy to reproduce experimental constraints, the configura- 
tion of the nuclei is less significantly altered and similar "norm" trajectories 
are produced with negligible consequences for core bounce [103]. 

We therefore conclude that the tight feedback between electron fractions and 
free proton abundances defines a limiting "norm" trajectory for the delep- 
tonization during infall. It is determined by the mass formula of nuclei ac- 
cording to its implementation in the equation of state. The deleptonization 
up to neutrino trapping can occur faster than "norm" by the availability of 
additional channels for electron capture, but hardly slower. The free proton 
mass fraction is not a relevant independent input quantity of the progenitor 
structure, it evolves rather symptomatically: in models that proceed along 
the "norm" trajectory, it is higher; in models that deleptonize faster, it is 
smaller. For example the most recent simulations with sophisticated nuclear 
input physics predict a deleptonization that is faster than the "norm" because 
captures on nuclei provide the most relevant electron capture channel through- 
out core collapse [21,96]. As we will discuss in more detail in section 3.3, the 
entropy also becomes shghtly smaller. Both effects contribute to reduce the 
proton abundances when treated consistently with the hydrodynamical evo- 
lution so that electron capture on free protons becomes an almost inactive 
channel. 

The large electron capture rates at densities higher than 5 x 10^^ g/cm^ are 
not relevant anymore because the neutrino scattering opacities get so large 
that the time scale for neutrinos to leave the star becomes larger than the 
production time scale. The neutrino phase space fills up and further electron 
captures are blocked for neutrino emission below the neutrino Fermi energy. 
In this regime the state of matter is best described by the density, entropy, 
and lepton fraction, Y^. The lepton fraction becomes the independent vari- 
able that slowly changes on a diffusion time scale while the electron fraction, 
"Ke; almost immediately represents the local equilibrium state between the 
fluid and the trapped neutrino radiation held. The equilibrium is established 
by neutrino-electron scattering. Due to the small electron mass, the energy 
transfer between neutrino and the relativistic electrons is always inelastic and 
efficient while the much more frequent scattering between the neutrinos and 
the more massive nucleons and nuclei only allows for momentum transfer — at 
least, of course, if we neglect for a moment the possibility of intrinsic excita- 
tions of nuclei in neutrino scattering events (see section 3.4 for a discussion of 
that aspect). 

Crucial for the total deleptonization during core collapse is the physics in the 
density range 5 x 10^° — 5 x 10^^ g/cm^ inbetween the two extremes, free escape 
and trapping. Figure 7 visualizes the interplay between neutrino production. 
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Fig. 7. Shown are the electron chemical potential, fig (straight solid line), and its 
distance to the difference between the nucleon chemical potentials, jj-e — IJ^n + IJ'p 
(dotted line). This quantity corresponds to the electron neutrino chemical poten- 
tial where the neutrinos are in equilibrium with matter (recognizable in the figure 
as the domain where the mean neutrino energy, E^, is equal to the equilibrium 
energy, EequH.)- The asterisks represent the average energy of neutrinos produced 
by electron capture. The neutrinos escape in a two-dimensional diffusion process 
in energy and space. Thick solid lines connect locations of optical depths r = 1, 
r = 10, and r = 100. For each emitted average energy, a thin solid line indicates 
the trajectory of fastest escape as a function of neutrino energy and density. At 
subnuclear densities (e.g. circle at p = 5 x 10^^ g/cm'^), low energy neutrinos start 
to escape and to free phase space for the downscattering of higher energy neutrinos 
around p = 10^^ — 10^^ g/cm^. In this region, the figure illustrates the nontrivial 
interplay between neutrino emission, thermalization, and diffusion. The left panel 
corresponds to a time slice at bounce in model G15 [31] with conventional input 
physics. The right panel corresponds to a simple test calculation where, in the spirit 
of Ref. [110], the orbit in the Bruenn scheme was unblocked to illustrate the 
discussion of the improved electron capture rates in section 3.3. 

neutrino thermalization, and neutrino diffusion in a snapshot at bounce in 
model G15 [31]. Figure 7 is built around the electron chemical potential, fie 
(thin solid line), which sets the energy scale for reactions with the degener- 
ate electrons. Also shown is the quantity /Xe — /^n + /^p (dotted line) where 
f^n — f^p is the difference between the neutron and proton chemical potentials. 
The quantity represents the neutrino chemical potential in regions where the 
neutrinos are in equilibrium with the fluid. Thick solid lines running from 
the upper left corner to the lower right indicate the optical depth, r, as a 
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function of density and neutrino energy. Contour lines for r = 1 (in aver- 
age one scattering before escape), r = 10 (in average 10 scatterings before 
escape, etc.), and r = 100 are shown. The contour hnes approximately run 
along E'^p^/^ — const, because of the increase of the scattering cross sections 
proportional to the squared neutrino energy [7]. The asterisks represent the 
average energy of emitted neutrinos by electron capture on free protons and 
on nuclei. For each asterisk in the graph, the trajectory of the most probable 
escape in the two-dimensional space of neutrino energy and matter density has 
been calculated and drawn with a thin sohd line. Note that the trajectories 
only connect the most probable of all possible displacements of test neutrinos 
in the density-energy space. They cannot represent a rate of actual displace- 
ments. The trajectories immediately confirm that neutrinos on the left hand 
side of the r = 1 line escape without further energy change while the neutrinos 
on the right hand side of the r = 100 line tend to scatter to the neutrino Fermi 
energy on a much smaller time scale than for the diffusion in space. The holes 
they leave in the equilibrium distribution will immediately be refilled by new 
neutrino emissions or scatterings. Note also the trajectory from the additional 
5 MeV test particle at 2 x 10^^ g/cm^, which can only scatter to the Fermi 
surface because the neutrinos are degenerate (The small steps in the Fermi 
energy are of numerical nature and due to the finite resolution of the energy 
phase space in the multi-group neutrino treatment). The other two 5 MeV 
test particles at 5 x 10^^ and 10^^ g/cw? demonstrate that at these densities 
very low energy neutrinos start to escape due to the relatively low scattering 
cross sections at these energies [51]. This frees phase space that is replenished 
from high energy neutrinos that are now able to down-scatter as indicated by 
the trajectories in the interesting density range between 5 x 10^° and 5 x lO"*^^ 
g/cm^. High energy neutrinos emitted at 10^^ g/cw? finally escape with a 
lower energy than neutrinos emitted at a lower energy at 10^^ g/cm'^. Hence, 
the most significant deleptonization relies on a two-dimensional diffusion pro- 
cess in the space spanned by neutrino energy and matter density and takes 
place during the transition from the density where neutrinos stream freely to 
the density where they are trapped. 

The physics of the core collapse does not depend on the evolution of the 
electron fraction alone. Also the entropy plays an important role. In general 
the entropy is low during the collapse as most of the energy is rather stored 
in internal excitation energy of the nuclei than in thermal energy of the gas 
due to the large partition functions for nuclei [44]. However, the entropy is 
changed by weak interactions in the core and has important consequences for 
the evolution (e.g., an increase in entropy will increase the number of protons 
available for electron capture). The entropy change As for a deleptonization 
AFe is given in Ref. [2] and references therein, 

TAS = -AFe (/^e - /^n + /^p) + Ag, (16) 
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where Ag is the energy transfer between matter and neutrinos. If we are 
not too strict about where in space the energy transfer between matter and 
neutrinos takes place (for neutrino-electron scattering rather along one of the 
trajectories in Fig. 7 than in one and the same fluid element), one may simply 
set Aq — /^YgEl^^^"^^ where E'^^'^^p^ is the average energy of the neutrinos when 
they escape the star [7]. This leads to the following expression for the entropy 
change (16): 

TAs = -AY, (/.e -l^n + fip- Et'^^^n ■ (17) 

Figure 7 now illustrates, that up to a density of 2 x 10^^ g/cm^ we find 
^escape ^ _ _|_ ^jjg eutropy in this regime is decreasing. Therefore, 
also the entropy might slightly contribute to the convergence to a "norm" 
trajectory during the deleptonization: If a fluid clement in one model has a 
higher temperature than in another model, more free protons are available for 
electron capture and the fluid element cools more rapidly toward the temper- 
ature in the other model. We also note that in this regime electron capture 
on nuclei produces neutrinos with lower escape energies and therefore con- 
tributes less to an entropy decrease than electron capture on free protons 
[111]. Between p = 2 x 10^^ g/cm^ and p = 2 x 10^^ g/cm^. Fig. 7 indicates 
that El^^"^^ < He — Pn + l^p such that the entropy increases because a part 
of the locally emitted neutrino energy is returned to the fluid by the ther- 
malization of the neutrino. At even higher densities, only the thermalization 
down to the neutrino chemical potential occurs on a fast time scale such that 
fast electron fraction changes do no longer cause entropy changes because of 
= Pe — l^n + l^p- The energy difference between p^ and E^^^^^^ is only re- 
turned to the fluid on the much longer diffusion time scale relevant for changes 
of the lepton fraction AY^. These expectations are consistent with the actual 
entropy evolution at the center shown in Fig. 10. 

3.3 Electron capture rates 

Calculations of the reaction rate for electron capture in the collapsing core 
requires two components: the appropriate electron capture reaction rates and 
the knowledge of the nuclear composition. The coupling of electron capture 
rates to energy-dependent neutrino transport adds an additional requirement: 
information about the spectra of emitted neutrinos. These spectra can be 
parametrized using the prescription of Ref. [112]. During the collapse most of 
the collapsing matter survives in heavy nuclei as the entropy is rather low [44] . 
Yg decreases during the collapse due to electron capture, making the matter 
composition more neutron rich and hence favoring increasingly heavy nuclei. 
Unlike in stellar evolution and supernova nucleosynthesis simulations, where 
the nuclear composition is tracked in detail via a reaction network [1,113], 
the composition used in supernova simulations is calculated by the equation 
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of state, which assumes nuclear statistical equilibrium (NSE). Typically, the 
information about the nuclear composition provided by the equation of state 
is limited to the mass fractions of free neutrons and protons, a particles and 
the sum of all heavy nuclei, as well as the identity of an average heavy nucleus, 
calculated either in the liquid drop framework [89] or based on a relativistic 
mean field model [114,115]. It should be noted that the most abundant nucleus 
is not necessarily the nucleus which dominate electron capture during the infall 
phase. For the evaluation of reaction rates on nuclei, due to the dependence 
on nuclear structure effects, a single nucleus approximation is not sufficient. It 
must be replaced by an ensemble average. Figure 8 shows the typical nuclear 
abundances [116] at two different conditions taken from a collapse trajectory 
from Ref. [117]. Before discussing how the detailed abundance information can 
be implemented in collapse simulations we briefly review the foundation for 
the calculation of electron capture rates during the infall phase. 

Traditionally, in collapse simulations the treatment of electron capture on nu- 
clei is schematic and rather simplistic. The nuclear structure required to derive 
the capture rate is then described solely on the basis of an independent-particle 
model for iron-range nuclei, i.e., considering only Gamow- Teller transitions 
from /7/2 protons to neutrons [44,2,119,3]. In particular, this model pre- 
dicts that electron capture vanishes for nuclei with neutron number N > 40, 
arguing that Gamow- Teller transitions are blocked due to the Pauli princi- 
ple, as all possible final neutron orbitals are already occupied in nuclei with 
> 40 [106]. These nuclei dominate the composition for densities larger than 
a few 10^° g cm^. As a consequence of the model apphed in previous collapse 
simulations, electron capture on nuclei ceases at these densities and the cap- 
ture is entirely due to free protons. It has been pointed out [107] that this 
picture is too simple and that the blocking of the Gamow- Teller transitions 
will be overcome by thermal excitations which either moves protons into the 
gg/2 orbit or removes neutrons from the p/-shell, in both ways unblocking the 
GT transitions. According to the work of Ref. [107], due to "thermal unblock- 
ing" GT transitions dominate again for temperatures of the order of 1.5 MeV. 
A more important unblocking effect, which is already relevant at lower tem- 
peratures is expected from the residual interaction which will mix the gQ/2 and 
higher orbitals with the pf shell [120,21]. 

The calculation of electron capture on nuclei during the collapse phase requires 
a model that is able to describe the correlations and at the same time the 
high density of levels that can be thermally populated at moderate excitation 
energies. Direct shell-model diagonalizations are not yet possible due to the 
large model spaces involved. The calculations can be done using the Shell 
Model Monte Carlo approach [121] which allows for the calculation of nuclear 
properties at finite temperature in unprecedentedly large model spaces. This 
model complemented with Random Phase Approximation calculations for the 
computation of the transitions necessary for the determination of the electron 
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Fig. 8. Abundances of nuclei for two set of conditions during the collapse. The upper 
panel represents typical conditions during the early collapse while the lower panel 
shows conditions near to neutrino trapping. A Saha-like NSE code has been used 
in the calculation of the abundances [118] (Adapted from [116]). 



capture rate has been used recently for the calculation of the relevant rates 
for nuclei in the mass range A = 65-112 [21]. Currently electron capture rates 
are available for all the nuclei in the mass range A = 17-112 and for the 
temperatures and densities relevant for the collapse [47,74,75,21,122]. 



Figure 9 compares the electron capture rates for free protons and selected 
nuclei along a stellar trajectory taken from [28]. These nuclei are abundant at 
different stages of the collapse. For all the nuclei, the rates are dominated by 
GT transitions at low densities, while forbidden transitions contribute sizably 
for p > 10^^ g cm~^. The electron chemical potential /ig and the reaction Q 
value are the two important energy scales of the capture process. (They are 
shown on the right panel of figure 9.) The value dependence of the capture 



21 



10= 
10^ 
10^ 

10° 
10-^ 
10-2 
10-3 
10-* 































'''''' 


---- 69^11 : 




..... 76Qai 




— - ^^Gel 


" 1 


89b|.^ 

1 ' 



10" 



10" 

p (g cm-^) 



10'= 




p (g cm ■>) 



Fig. 9. (left panel) Comparison of the electron capture rates on free protons and 
selected nuclei as function of density chemical potential along a stellar collapse 
trajectory taken from [28]. (right panel) energy scales relevant for the determination 
of the electron capture rates, fie is the electron chemical potential. (Q) and Qp are 
the average Q value for electron capture on nuclei and protons respectively. 



rate can be well approximated by [48]: 
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where x = Q/T, r] = (/ig — Q)/T, K = 6146 s, and B represents a typical 
(Gamow- Teller plus forbidden) matrix element. The quantities Fk are the rel- 
ativistic Fermi integrals of order . Q represents the effective Q- value for the 
capture rate including thermal effects. In equation (18), we use the convention 
that Q is positive for capture on neutron rich nuclei and protons. For low den- 
sities (< 10^° g cm~'^), fie ^ Q and the term F2 dominates. The rate is then 
larger for the nuclei with smaller Q values. For intermediate densities (10^°- 
10^^ g cm"'^), the term F3 dominates and the rate still has some dependence 
on the Q- value. For high densities (> 10^^ g cm~^), fie^ Q so that the term 
F4 dominates and the rate becomes independent of the Q-value, depending 
only on the total GT strength, but not its detailed distribution. 

According to figure 9, the electron capture rate on a proton is larger than that 
for individual nuclei. However, this is misleading as the low entropy keeps the 
protons significantly less abundant than heavy nuclei during the collapse. Sim- 
ulations of core collapse require then the knowledge of the detailed abundances 
of all the nuclei present. As the commonly used equations of state [89,114,115] 
do not provide this detailed information, a Saha-like NSE was used for the 
calculation of the abundances in refs. [21,96]. Once the abundances are con- 
sidered the reaction rate for electron capture on heavy nuclei {Rh = J2i Y^iK, 
where the sum runs over all the nuclei present and Yi denotes the number 
abundance of species i) dominates over the one of protons {Rp = YpXp) by 
roughly an order of magnitude throughout the collapse [21,96]. The first sim- 



The fermi integrals are defined by: -Ffe(^) = /o°° dxx^ / (1 + exp(x — 77)) 
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ulations that considered electron capture on nuclei with N > 40 were probably 
the ones of Rcf. [111]. These authors determined the electron-capture rate by 
detailed balance from their computed absorption rate. This corresponds to 
electron capture on the thermally populated Fermi and GT_ resonances in 
the parent nucleus shown in figure 1 (also called backresonances according to 
the notation of [45]) that decay to the daughter ground state. This backreso- 
nance contribution is in general, depending on the nucleus, a factor between 
10 and 100 smaller than the direct transition from the parent ground state to 
the daughter Gamow- Teller resonance. This explains why, in contrary to the 
results of refs. [21,96], electron capture on nuclei was too small in Ref. [Ill] so 
that it was concluded that electron capture on protons dominates over capture 
on nuclei. 




pc (g cm ^) Pc (g cm ) pc (g cm ) 

Fig. 10. Comparison of the evolution of several quantities at the center of a 15 Mq: 
Yg is the number of electrons per baryon, Y[ep is the number of leptons per baryon s is 
the entropy per baryon, T is the temperature, dE/dt is the neutrino energy emission 
rate per baryon, {Ei,) is the average energy of the emitted neutrinos. An abrupt 
change in the emission occurs when the domain where the LMSH rates are specified 
is exceeded at a density of 10^^ g/cm^. The initial presupernova model was taken 
from [76] . The thin line is a simulation using the Bruenn parametrization [2] while 
the thick line uses the LMSH rate set (sec text). Both models include the general 
relativistic effects discussed in section 3.1 and were calculated by the Garching 
collaboration (Courtesy of M. Rampp and H.-Th. Janka). 

The effects of the more realistic implementation of electron capture on heavy 
nuclei have been evaluated as described above in independent self-consistent 
neutrino radiation hydrodynamics simulations by the Oak Ridge and Garching 
collaborations [96,123]. The results obtained by the Garching collaboration are 
shown in Fig. 10 [123] for the evolution of several quantities at the center of 
a 15 M0 star using the standard treatment of Bruenn [2] and the new rates 
for heavy nuclei (denoted LMSH). With the improved treatment of electron 
capture rate on heavy nuclei the total electron capture rate (heavy nuclei 
plus protons) is larger than in the Bruenn treatment resulting in smaller Yg 



23 



and liepton values. The left panel shows the evolution of these quantities at 
the center of the star as functions of the increasing density during collapse. 
The small oscillations in the trajectories above 5 x 10^^ g/cm^ are due to the 
same numerical effect of finite energy resolution under neutrino-degenerate 
conditions that we have already seen in Fig. 7. The panel demonstrates that 
the Bruenn treatment leads to a faster deleptonization for representative nuclei 
with < 40 (at low density) while this channel is blocked afterwards such that 
the further deleptonization is delayed and falls short of the LMSH treatment 
where electron captures on nuclei continue. The panel on the right shows the 
energy per baryon and second emitted in neutrinos at the center of the star. 
The inset shows the mean energy of the produced neutrinos. With the LMSH 
treatment, more neutrinos are produced with lower mean energy than with 
the Bruenn treatment. The center panel shows the entropy and temperature 
evolution at the center of the star. The general entropy trajectory is consistent 
with the discussion following Eq. (17). The entropy is slowly descending at low 
density because _E'^scapc -g [g^Ygci than fie — fin + fJ'p- As soon as neutrino-electron 
scattering becomes effective to downscatter neutrinos, i^^scapc ^^^[s below above 
chemical potential difference (see Fig. 7, left panel) and the transferred heat 
starts to increase the fluid entropy until the diffusion time scale becomes larger 
than the dynamical time scale. With the LMSH treatment one would now 
be tempted to expect that the central entropy would rise to larger values 
than with the Bruenn treatment because the deleptonization AYe in Eq. (17) 
increases according to the values shown in the left panel. Additionally the 
neutrinos are emitted with lower energies so that £'^^='^p^ might also decrease 
somewhat to favor a steep entropy increase. However, the center panel shows 
clearly that the LMSH treatment leads to smaller central entropies. In order 
to better illustrate speculations about the possible reasons of this interesting 
result, we performed a simple test calculation for a comparison in Fig. 7 where 
one state in the neutron orbit was kept unblocked in the Bruenn treatment. 
One immediately sees in the right panel that the average energy of emitted 
neutrinos (asterisks) is lower than with the Bruenn treatment and that less 
neutrino-electron scattering takes place before the escape. This difference in 
the thermalization path reduces the diffusion time. For an understanding of 
the entropy differences however, it is crucial to realize that the difference in 
the deleptonization between the Bruenn and LMSH rates occurs before the 
density is reached where the entropy rises. We define a first phase, reaching 
from low densities to p = 3 x 10^'^ g/cm^, so that the net deleptonization 
for the two rates is exactly the same (left panel in Fig. 10). As both rates 
include electron capture on nuclei at these densities, Eq. (17) does not predict 
noticeable differences in the entropy evolution. In a second phase, reaching 
from p = 3 X 10^° to p = 2 x 10^^ g/cm^, electron capture on nuclei is blocked 
in the Bruenn treatment and a significantly smaller net deleptonization —AYe 
in Eq. (17) occurs than with the LMSH rates. However, Fig. 7 shows that in 
this phase the second factor, /ig — /in + iJ^p — E^^^^^"^^, is much smaller than in 
the LMSH case because electron capture on free protons produces neutrinos 
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with larger escape energies £;^scape ^^j-^g^j^ electron capture on nuclei. The two 
effects cancel in the evaluation of TAs in Eq. (17) and the entropy evolution 
in this phase is again similar with both sets of rates. Now, in the third phase, 
reaching from p = 2 x 10^^ to p = 2 x 10^^ g/cm^, there is no significant 
difference in the net deleptonization for the two cases anymore, —AYe ~ 0.09, 
but the lower Ye values increase the difference between the neutron and proton 
chemical potentials, — //p, in the LMSH case. On the other hand, the energy 
of escaping neutrinos, £;^scape^ decreases, but only about half as much. Both 
effects together cause a difference ~ 2 MeV in the factor /Xg — /^n + /^p — 
^escape ^ which seems consistent with the displayed difference in the entropy 
evolution before trapping in the center panel of Fig. 10. The increased neutron 
to proton asymmetry in the LMSH case also leads to heavier nuclei with higher 
opacities for a given density and neutrino energy (see section 3.5). This shifts 
the deleptonization region to slightly lower densities where the energy scales 
are smaller. This might also contribute to the entropy difference. 



3.4 Inelastic Neutrino-Nucleus Interactions 

Although the classical thermalization mechanism for the neutrinos proceeds 
through inelastic scattering on electrons, Eq. (8), neutrino-induced reactions 
on nuclei can also occur [124,111]. Charged-current {ugjC') reactions on nu- 
clei are unimportant during the collapse compared with inelastic neutrino- 
electron scattering due to the strong Pauli blocking of the final electron phase 
space [111]. This is so even if finite temperature effects are considered [125]. 

However, according to Rcf. [Ill], neutrino- nucleus inelastic scattering plays 
an as important role as neutrino-electron scattering. The infiuence of finite 
temperature was studied in Ref. [126] using an independent particle model. 
More recently, neutrino cross sections for finite temperature have been calcu- 
lated using shell-model GT distributions [127,128]. This approach has been 
validated for the calculation of neutrino-nucleus inelastic scattering against 
high resolution magnetic dipole data [129]. Finite temperature enhances the 
cross section for low energy neutrinos. This is due to the fact that at finite 
temperatures the initial nucleus can reside in excited states which can be 
connected to the ground state by large GT matrix elements. These states are 
deexcited by inelastic neutrino scattering. In this case, contrarily to what hap- 
pens in neutrino-electron scattering, the final neutrino energy is larger than 
the initial. Currently no hydrodynamical collapse simulations are available 
that include both electron capture on nuclei and neutrino-nucleus inelastic 
scattering. Nevertheless, the simple fact that the cross sections are of sim- 
ilar magnitude than scattering on electrons [127] suggest their inclusion in 
supernova neutrino transport codes. 
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Other neutrino processes such as nuclear deexcitation by neutrino pair pro- 
duction, Eq. (11), have been discussed in Ref. [126], but the estimated rates 
are probably too small for these processes to be important during the collapse. 



3.5 Neutrino opacities 

After the detailed description of the new electron capture rates in section 3.3 
and the emphasis of the importance of neutrino thermalization in sections 3.2 
and 3.4, we conclude the discussion of nuclear input for core-collapse models 
with a few pointers to the third relevant ingredient: neutrino opacities. The 
fundamental neutrino opacity in core collapse is provided by neutrino scatter- 
ing on nucleons. Depending on the distribution of the nucleons in space and 
the wavelength of the neutrinos, various important coherence effects can oc- 
cur: Most important during collapse is the binding of nucleons in nuclei with 
a density contrast of several orders of magnitude to the surrounding nucleon 
gas. For elastic neutrino-nucleus scattering one usually makes the simplifying 
assumption that the nucleus has a J = 0"*" spin/parity assignment. This is 
appropriate for the ground state of even-even nuclei. The scattering process 
is then restricted to the Fermi part of the neutral current (pure vector cou- 
pling) [130,131]. Because of coherent scattering, the cross section scales with 
A^, except from a correction ~ (A^ — Z)/A arising from the neutron excess. 
This assumption is, in principle, not correct for the ground states of odd-A 
and odd-odd nuclei and for all nuclei at finite temperature, as then J > 
and the cross section will also have an axial-vector Gamow-Teller contribu- 
tion. However, the relevant GTq strength is not concentrated in one state, but 
rather fragmented over many nuclear levels. Thus, one can expect that the 
GT contributions to the elastic neutrino-nucleus cross sections are in general 
small enough to be neglected. A useful comparison of inverse mean free paths 
at the important density p — 10^^ g/cm^ is given in [132]. It is found that 
K+n/K+A ~ 3 X 10-2, Ai.+e/A^+A ~ 2.5 X 10"^ X^+He/K+A ~ 10~^, and 
X^+p/X^+A ~ 5 X 10"^ 

Further corrections are necessary: with an increasing ratio between the Coulomb 
potential of the positively charged ions and their thermal energy, the average 
separation between nuclei will more strongly peak around the value of most 
efficient packing. The neutrino opacities are then to be corrected by an ion- 
ion correlation function {Si^n {Ey)) < 1 [133,134,135]. Its consideration in core 
collapse simulations lowers the trapped lepton fraction at bounce by 0.015 
and increases the central entropy per baryon by 0.12 kB [132]. It is found that 
the sizable entropy increase is not only due to the increased deleptonization, 
but also to the fact that the correlation effect is most pronounced for the low 
energy neutrinos with a long wavelength. Hence the entropy increase is addi- 
tionally supported by a decrease of E^^^^ in the discussion following Eq. (17). 
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As current core collapse models proceed toward including the full ensemble of 
nuclei instead of just a representative nucleus (cf. Fig. 8), it becomes rather 
non-trivial how to adequately determine correlation effects in the ion mixture, 
see for example [134,136]. 

The situation is even more interesting in the phase transition from isolated 
nuclei to bulk nuclear matter where the nuclei or the holes inbetween them 
are strongly deformed. Various pasta-like shapes may be assumed. Correla- 
tion effects in this phase could also affect the neutrino opacities [99]. For an 
immediate effect on core collapse, however, it would be required that this 
"pasta-phase" would reach to fairly low densities in order to affect the opac- 
ities at the neutrinospheres where the neutrino luminosities and spectra are 
set (cf. Fig. 7). 

A more detailed description of the neutrino opacities in nuclear matter is given 
elsewhere in this volume [6] . An extensive quantitative overview of the rates of 
most reactions has been provided in Ref. [111]. Useful fitformulae are collected 
in Ref. [134]. 



4 Effects on the postbounce evolution 

Stellar core collapse and supernova explosion have often been discussed as 
one and the same event. This view has also observationally been confirmed 
by the detection of few neutrinos from supernova 1987A [137,138,139]. The 
connection has further been nourished by the early theory of a prompt ex- 
plosion mechanism where, as a consequence of bounce, a strong shock rushes 
through the outer core to eject the outer layers. However, we are not aware 
of any simulation with reasonable neutrino physics that would have predicted 
a prompt explosion since the Lattimer-Swesty equation of state [89] has be- 
come standard. In contrary, due to the stiff behaviour of the equation of state 
around nuclear density the initial shock is weak as discussed in section 3.1. 
Additionally it starts at a very deep mass coordinate because of the efficient 
deleptonization discussed in sections 3.2-3.3. In fact, recent spherically sym- 
metric simulations of the postbounce evolution predict that the shock stalls 
already ~ 3 ms after bounce, i.e. before or at the time the neutrino burst 
is launched from neutrino-transparent layers and ~ 40 ms before a gain layer 
can develop for neutrino heating to become effective. Many publications of su- 
pernova models include figures with shock trajectories that show a continued 
expansion of the shock front to about 150 km at 100 ms after bounce. Note, 
however, that there is no valid concept of shock energy that reaches beyond 
the few milliseconds after which the shock has stalled. The continued expan- 
sion of the accretion front over at least 90% of the ~ 100 ms interval is due 
to the accumulation of material on the protoneutron star: heavy nuclei from 
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the outer layers are accelerated in the gravitational potential until they fall 
upon the accretion front where their kinetic energy is converted into heat. The 
accretion front is displaced to a larger radius mainly by the increased volume 
of the accumulated hot and dissociated matter and not by shock energy. 

The sophisticated nuclear physics required in core collapse models and the 
complicated dynamics in the ensuing supernova explosion suggest a point of 
view where the theory of core collapse may well be as detached from supernova 
explosions as, for example, the dynamics of star formation from stellar evolu- 
tion. On the one hand, core collapse involves the nuclear and weak interaction 
physics of heavy nuclei under electron degenerate conditions in layers that are 
unstable to gravitational collapse. The supernova explosion, on the other hand, 
involves fundamentally multi-dimensional neutrino-radiation-hydrodynamics 
in hot dissociated matter that is in convective motion around a possibly mag- 
netized compact object. Also the technical refinement of simulations of the 
collapse phase has always been roughly 10 years ahead of the corresponding 
simulations of the postbounce phase. Based on the current understanding, 
it would be an optimistic expectation that any very relevant change in one 
domain (core collapse) would automatically induce a revolution in the other 
(supernova explosion). 

Important changes have been discovered in the core collapse event. In con- 
ventional simulations, electron capture on nuclei has been treated as blocked 
for nuclei with iV > 40. Recent models that include realistic electron capture 
rates beyond this limit show that electron capture on these nuclei are not 
only relevant, but even provide the dominant channel for the deleptonization 
during collapse [21,96], with important consequences for core-bounce: In the 
interior of the protoneutron star the electron fraction is reduced by ~ 10% 
when compared to models that don't include electron captures on nuclei. This 
change corresponds to a reduction of 20% in the mass of the homologous core 
and translates directly to a reduction of the mass interior to the point of shock 
formation. Moreover, the shock is launched with a ~ 15% smaller velocity dif- 
ference. The weaker shock with more iron core overlying it turns even earlier 
(and at a deeper mass coordinate) into an expanding accretion front than in 
models with "standard" physics. Furthermore, the improved treatment of elec- 
tron capture also affects the outer layers where the average neutron number 
is well below = 40. For these conditions the treatment of electron cap- 
ture by Bruenn [2] results in more electron capture than with the shell-model 
based electron capture rates [20] . The smaller neutronization in the outer lay- 
ers slows the collapse in these layers which further diminishes the growth of 
the electron capture rate. These changes delay the accretion of matter and 
reduce the ram pressure opposing the expansion of the accretion front. Both 
effects can modify the trajectory of the accretion front as a function of time. 

In spherically symmetric simulations the accretion front reaches a maximum 
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expansion ~ 100 ms after bounce when the compression of matter in the 
coohng region starts to exceed the volume of accumulated shock- and neutrino- 
heated matter at the accretion front. However, because the shock stalls within 
milliseconds in old and new models, the significant changes in the core collapse 
physics cannot dynamically push through to the explosion phase. Also the 
temporal variations in the accretion rate tend to average over ~ 100 ms, so 
that the maximum radius in the postbounce expansion of the accretion front is 
not significantly sensitive to the improvements of the electron capture rates on 
nuclei. Nevertheless, core collapse still sets the stage for a delayed supernova 
explosion and several consequences of the nuclear physics during core collapse 
survive for a longer time than the often overemphasized shock dynamics, or 
are in principle directly observable. 

Due to the extended accretion phase of delayed explosions, most if not all 
of the iron core is packed in the protoneutron star. Hence, differences in the 
physics of core collapse will most likely translate to differences in the interior 
of the protoneutron star. The deleptonization during collapse and the early 
shock energetics leave a clear imprint on the entropy profile at high densi- 
ties. Indeed, models that include electron capture on nuclei show considerably 
different entropy and lepton fraction gradients when compared to standard 
treatments of electron capture [96]. These gradients are responsible for the 
location, extent and strength of the protoneutron star convection and other 
possible fluid instabilities. Seminal simulations by [140] favored a strong cou- 
pling between the protoneutron star interior and the exterior layers subject 
to neutrino heating. The equation of state used at that time suggested fluid 
instabilities which, implemented by a mixing-length approach in spherically 
symmetric simulations, supported the existence of delayed neutrino- driven su- 
pernova explosions. Analytical investigations based on the Lattimer-Swesty 
equation of state [141] and recent numerical simulations [36] find convectively 
stable neutrinospheres that efficiently block a rapid release of neutrinos from 
the convectively unstable protoneutron star interior. But a change in the nu- 
clear input physics always carries the potential to affect the structure of the 
protoneutron star in a way that possibly revives the strong coupling of its in- 
terior to the outer layers. A recent detailed investigation of the protoneutron 
star stabihty is given in Ref. [37]. 

The neutrino luminosities and energies during collapse, bounce, and the inter- 
esting time of the launch of the electron neutrino burst are in principle more 
directly observable. Figure 11 shows the luminosity and root-mean-square en- 
ergy of the emitted electron neutrinos and antineutrinos (sampled at 500 km 
radius in the comoving coordinate frame) between 100 ms before bounce and 
100 ms after bounce. One model (thick lines) includes electron captures on 
nuclei (LMSH) and the other (thin lines) uses the standard Bruenn treat- 
ment. The luminosity shows a clear delay (2 ms) in the breakout burst due to 
a deeper launch of the shock in the LMSH case. Consistent with the depen- 
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Fig. 11. The neutrino luminosity and root-mean-square energy (at 500 km) as a 
function of time from bounce for a 15 Mq model. The thin lines show this evolution 
for a simulation using the Bruenn parametrization [2] , while the thick lines show this 
evolution for a simulation using the LMSH rates [21,96]. The solid lines correspond 
to electron type neutrinos; the dashed lines correspond to electron type antincutri- 
nos. Both cases were calculated in a general relativistic model by the Oak-Ridge 
collaboration (adapted from [96]). 



dence of the neutrino burst intensity on the shock strength [31], the LMSH 
neutrino burst is also more extended in time. Before bounce the neutrino lu- 
minosities are similar but the average energies arc smaller in the LMSH case 
due to the larger Q-values involved for electron capture on nuclei as discussed 
above. After bounce, during the first 50 ms the LMSH model emits ~ 15% 
more energy than the Bruenn model. This is mainly the result of a larger mean 
electron neutrino energy in the LMSH model due to a deeper and hotter neu- 
trinosphere [96]. 

Neutrino-induced reactions on nuclei are expected to contribute only modestly 
to the shock revival, due to the small abundance of nuclei during this phase. 
It has been suggested that the shock revival can be aided by "preheating" of 
matter ahead of the shock by absorption of electron neutrinos emitted in the 
breakout burst [124]. These neutrinos can partly dissociate the matter, mainly 
iron and sihcon isotopes, before the arrival of the shock. The early explorations 
of Ref. [Ill] found no significant preheating of unshocked matter. Once the 
matter is processed by the shock, large amounts of ^He are present just below 
the stagnated shock. The possibility of reheating this material by inelastic u^- 
^He scattering, where i^^ denotes muon or tau neutrino and antineutrino, or 
heating the material ahead of the shock by i/a;-nucleus inelastic scattering was 
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also investigated in Ref. [Ill] with no major effect found. Nevertheless, only 
recently reliable neutrino- nucleus cross sections have become available [128] 
that allow for a complete study of the influence of neutrino-nucleus scattering 
both in the collapse and postbounce evolution. 



5 Conclusions 

Core-collapse supernovae present a rich interplay between nuclear physics and 
neutrino radiation (magneto-)hydrodynamics. Nuclear physics sets the scene 
in which different hydrodynamical instabilities could develop that eventually 
may produce robust explosions by the neutrino-driven mechanism (see for 
example refs. [142,143,39] for recent overviews over the activities of collab- 
orations engaging in accurate neutrino transport). The nuclear equation of 
state governs stellar core collapse, core bounce and neutron star formation. 
The complicated interplay between deleptonization, neutrino thermalization, 
and diffusion based on different weak interaction processes — the most impor- 
tant of them is electron capture on nuclei — determine the collapse dynamics, 
the position of shock formation and the structure of the layers into which 
the shock expands. Advances in the calculations of these processes coupled 
to improvements in one and multi-dimensional accurate neutrino transport 
simulations has allowed to gain further insight to the physics of the collaps- 
ing core. Future advances on the nuclear physics side should be guided by 
the advent of new radiative-ion beam facilities that will help to constrain the 
nuclear models used for the determination of the equation of state and the dif- 
ferent nuclear weak-interaction processes. These advances coupled to current 
efforts to develop multidimensional magneto-hydrodynamical codes with accu- 
rate neutrino transport will foster a definitive and quantitative understanding 
of supernova explosions. 
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